특성 다항식 $p(\lambda) = \det(A - \lambda I)$는 고유값을 정의하는 이론적 기초이지만, 고차원 시스템에서는 수치적으로 '불안정'하고 계산적으로 비효율적입니다. 실제 응용 사례(예: 파동 전파를 위한 슈튀름-리우빌 시스템 해결)에서 다항식 근의 계수 변화에 대한 민감성 때문에 직접 전개는 보조적인 선택입니다.
연속파에서 이산 행렬로
현수 또는 막의 진동은 파동 방정식으로 결정됩니다:
$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$
해 $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$를 찾기 위해 우리는 다음을 해결해야 합니다 슈튀름-리우빌 시스템:
$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$
이산화의 복잡성
연산자의 이산화는 $Aw = -0.04 \frac{\rho}{p} \lambda w$와 같은 행렬 방정식을 유도합니다. $4 \times 4$ 삼각행렬의 경우 $p(\lambda)$는 다룰 수 있지만, 메시가 세밀해지면서($n$ 증가), 우리는 두 가지 장벽에 부딪힙니다:
- 아벨-루피니 제한: 다항식의 근에 대해 $n \ge 5$인 경우 대수적 해가 존재하지 않습니다.
- 반올림 민감성: 고차원 시스템에서 하나의 항목이 $10^{-10}$ 소수 자리에서 변경되면 고유값이 주문 단위로 바뀔 수 있습니다(윌킨슨 다항식 현상).
수치적 필요성과 전문 라이브러리
전문 수치 라이브러리(имсл, НАГ)는 원시 특성 다항식을 피하고 대신 근사치를 구하기 위한 반복적 절차를 사용합니다:
- IMSL 라이브러리: 선형 최소제곱, 세제곱 스플라인, 그리고 빠른 푸리에 변환을 사용합니다.
- NAG 라이브러리: 최소제곱 다항식 근사 및 $l_1/l_{\infty}$ 기준 적합을 사용합니다.
시스템 $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$의 고유값을 근사할 때, 근 찾기보다는 이산 최소제곱과 반복적 탐색에 의존합니다.
🎯 이론적 도구와 수치적 위험
특성 다항식 $p(\lambda) = \det(A - \lambda I)$는 증명에 필수적이지만 계산에는 위험합니다. 물리학의 실제 고유값 문제는 안정성을 유지하는 반복 변환(예: QR)을 통해 해결됩니다.